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ABSTRACT 

A new dual-frequency (Ku and Ka band) nadir-pointing Doppler radar on the high-altitude NASA ER-2 
aircraft, called the High-Altitude Imaging Wind and Rain Airborne Profiler (HIWRAP), has collected data 
over severe thunderstorms in Oklahoma and Kansas during the Midlatitude Continental Convective Clouds 
Experiment (MC3E). The overarching motivation for this study is to understand the behavior of the dual- 
wavelength airborne radar measurements in a global variety of thunderstorms and how these may relate to 
future spaceborne-radar measurements. HIWRAP is operated at frequencies that are similar to those of the 
precipitation radar on the Tropical Rainfall Measuring Mission (Ku band) and the upcoming Global Pre- 
cipitation Measurement mission satellite’s dual-frequency (Ku and Ka bands) precipitation radar. The air- 
craft measurements of strong hailstorms have been combined with ground-based polarimetric measurements 
to obtain a better understanding of the response of the Ku- and Ka-band radar to the vertical distribution of 
the hydrometeors, including hail. Data from two flight lines on 24 May 2011 are presented. Doppler velocities 
were —39 m s _1 at 10.7-km altitude from the first flight line early on 24 May, and the lower value of —25 m s _1 
on a second flight line later in the day. Vertical motions estimated using a fall speed estimate for large graupel 
and hail suggested that the first storm had an updraft that possibly exceeded 60 m s~ 1 for the more intense part 
of the storm. This large updraft speed along with reports of 5-cm hail at the surface, reflectivities reaching 
70 dBZ at S band in the storm cores, and hail signals from polarimetric data provide a highly challenging 
situation for spaceborne-radar measurements in intense convective systems. The Ku- and Ka-band re- 
flectivities rarely exceed —47 and —37 dBZ, respectively, in these storms. 


1. Introduction 

Currently flying spaceborne precipitation and cloud 
radars, including the Ku-band precipitation radar (PR) 
on the Tropical Rainfall Measuring Mission (TRMM) 
(Kummerow et al. 2000) and the W-band CloudSat 
(Stephens et al. 2008), provide global coverage of a 
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variety of weather and cloud systems. With the success 
of these radars, the next spaceborne meteorologically 
oriented radars to be launched are the Global Pre- 
cipitation Mission (GPM) dual-frequency (Ku and Ka 
band) precipitation radar (DPR) in 2014 (Hou et al. 
2008) and the W-band Doppler EarthCare Cloud Pro- 
filing Radar (Horie et al. 2010) in 2016. Because of 
spaceborne-antenna size constraints, the current and 
future radars all use relatively short wavelengths that 
pose challenges for observing deep convective storms 
with high rain rates or large ice particles including 
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graupel and hail. As the hydrometeor sizes become 
comparable to or larger than the radar wavelength, Mie 
scattering occurs rather than Rayleigh scattering. As 
a consequence, the radar returns from typical ground- 
based S-band radars and those from shorter-wavelength 
satellite radars can be drastically different. TRMM and 
CloudSat missions have provided a wealth of informa- 
tion on intense deep convection, but we do not yet have 
comparable Ka-band measurements except through 
simulations (e.g., Tanelli et al. 2011). 

Several studies have documented global characteris- 
tics of extremely intense deep convection (e.g., Zipser 
et al. 2006) by using radar profile characteristics with 
echo height (such as the 50-dBZ contour) as a surrogate 
for updraft strength, which is not directly measured by 
TRMM. These analyses have provided the ability to 
quantify the intensity of convection globally and to un- 
derstand a host of extremely important properties of the 
convection such as their regional and diurnal variations. 
It is also well known that when radio waves propagate 
through regions of high liquid water content or through 
regions with large liquid or ice hydrometeors they can be 
strongly attenuated. In CloudSat measurements of deep 
convection, large portions of high-reflectivity cores are 
often entirely missing because of severe W-band atten- 
uation (e.g., Mitrescu et al. 2008). The strong attenua- 
tion in deep convection could also affect the TRMM PR 
at Ku band. It is often difficult to quantify hydrometeor 
characteristics in convection with single-frequency ra- 
dars such as TRMM or CloudSat. Furthermore, it is dif- 
ficult to know whether graupel or hail is present and the 
extent of Mie effects, attenuation, and even multiple 
scattering (e.g., Battaglia et al. 2006; Battaglia et al. 2010). 

There are several downward-looking airborne radars 
that simulate spaceborne cloud and precipitation mea- 
surements in deep convection, hurricanes, and other 
weather phenomena. The only radar measurements to 
capture the full depth of deep convection have been with 
the single-frequency (X band) ER-2 Doppler Radar 
(EDOP) on a National Aeronautics and Space Admin- 
istration (NASA) ER-2 aircraft (e.g., Heymsfield et al. 
2010). The dual-frequency (Ku and Ka band) Airborne 
Precipitation Radar (APR-2) on a medium-altitude 
(<11 km) DC-8 aircraft has provided measurements in 
hurricanes and occasional strong tropical oceanic con- 
vective cells (e.g., Sadowy et al. 2003; Durden et al. 2003), 
but it is limited primarily by considerations of aircraft 
safety when studying land-based deep convection. 

The motivation for this paper is to understand the 
behavior of dual-frequency Ku-/Ka-band airborne radar 
measurements in deep convection through case study 
and how these measurements may relate to spaceborne- 
radar measurements. One question that needs to be 


addressed for the upcoming GPM DPR is how well the 
Ku-/Ka-band radar performs in strong midlatitude 
convection. The TRMM PR has observed severe con- 
vective events within ±35° latitude, such as in the 
southern portion of the U.S. southern Great Plains 
(SGP). There have been no dual-wavelength Ku-/Ka- 
band airborne or spaceborne measurements from these 
more intense storms, however. The Ka-band wavelength 
on DPR is better in more moderate rain and snow since 
it suffers from attenuation in convection. In this paper, 
we will present the first dual-frequency Ku- and Ka- 
band radar measurements over severe hailstorms from 
the High-Altitude Imaging Wind and Rain Airborne 
Profiler (HI WRAP). The measurements were taken 
near central Oklahoma during the Midlatitude Conti- 
nental Convective Clouds Experiment (MC3E). We 
begin with a HIWRAP instrument description and data 
processing in section 2, followed by an overview of the 
two severe hailstorms and their environmental condi- 
tions in section 3. Section 4 presents HIWRAP and 
ground-based polarimetric observations of these hail- 
storms that have significant scattering and attenuation at 
Ku and Ka bands. Section 5 presents interpretations 
from the observations and their implications for TRMM, 
GPM, and other future satellite missions. To be specific, 
the magnitudes of updrafts and the usefulness of the dual- 
frequency ratio (DFR) are examined. Since knowledge of 
hydrometeor size distributions is poor for the strong 
convective regions studied, simple single-hailstone cal- 
culations are presented to provide insight on DFR. 

2. HIWRAP instrument description and processing 

HIWRAP is a dual-frequency (Ku and Ka band) and 
dual-beam conically scanning (30° and 40° incidence an- 
gles) Doppler radar system that was developed for remote 
operation on the NASA Global Hawk for studying winds 
in precipitation regions in hurricanes (Li et al. 2011). It 
utilizes solid-state power amplifiers, versatile waveforms, 
and pulse compression to achieve simultaneous transmit 
and receive at different frequencies and different beam 
positions. It flew for the first time on the NASA Global 
Hawk during the Genesis and Rapid Intensification Pro- 
cesses field campaign in 2010 with a downward-looking 
conical scan antenna (Braun et al. 2013). HIWRAP has 
frequencies that are similar to those of GPM, and it flew 
during MC3E with a new nonscanning (nadir pointing) 
antenna designed for the ER-2 (Table l). 1 For each radar 


1 The 33.7-GHz frequency was used for MC3E rather than 
HIWRAP’s upper band (35.5 GHz), which is similar to the GPM 
frequency, because the ER-2 antenna had better voltage standing 
wave ratio performance at the lower sideband. 
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Table 1. HIWRAP ER-2 system parameters during MC3E. 


Parameters 

Ku band 

Specifications 

Ka band 

Frequency (GFtz) 

Chirp mode: 13.915 
Pulse mode: 13.904 

Chirp mode: 33.733 
Pulse mode: 33.716 

Transmitter peak power (W) 

25 

8 

Antenna gain (dB) 

34.8 

42.0 

Antenna 3-dB beamwidth (along track; °) 

3.07 

1.23 

Antenna 3-dB beamwidth (cross track; °) 

2.96 

1.19 

PRF (Hz) 

4516/3589 

4516/3589 

Pulse width (/rs) 

Chirp mode: 20 
Pulse mode: 2 

Chirp mode: 20 
Pulse mode: 2 

Receiver bandwidth (MHz) 

2 

2 

Doppler range (ms -1 ) 

±97 

±40 

Min detectable reflectivity (dBZ) (at 10 km and 0.5-s avg) 

Chirp mode: —0.5 
Pulse mode: 9.5 

Chirp mode: -10.6 
Pulse mode: —2.9 


frequency (Ku or Ka band), a pulse sequence of 2, 20, 
and 2 /as was used for the data in this paper (McLinden 
et al. 2013). The first 2-/rs short pulse is used to obtain 
returns near the surface since strong surface returns 
contaminate the near-surface rain signal through pulse 
compression range sidelobes in the chirp channel. The 
second short pulse is used to obtain returns near the radar 
in the “blind” zone of the 20-/rs chirp pulse. HlWRAP’s 
beams are not matched as would be desired for dual- 
frequency algorithms since matched beams would 
decrease the Ka-band radar sensitivity significantly 
(~8 dBZ). Some of the effects from beam mismatch can 
be reduced through averaging the higher-resolution Ka- 
band measurements along track, but such averaging was 
not performed here. 

Calibration was performed by monitoring radar sys- 
tem stability with an internal calibration during flight 
and then by using the ocean-surface normalized radar 
cross section (NRCS) for external calibration (e.g., 
Tanelli et al. 2006). The radar system stability is moni- 
tored since it undergoes extreme temperature changes 
during the flight. This is accomplished using a calibration 
pulse obtained from an attenuated sample of the trans- 
mit pulse, which is passed through the digital receiver 
along with the received weather returns. An initial cal- 
ibration was performed using this sampled calibration 
pulse along with the measured variables in the radar 
equation. This provides a continuous end-to-end system 
calibration during the flight. The final calibration was 
performed using the ocean surface as an external ref- 
erence, similar to previous literature approaches. As 
mentioned in Tanelli et al. (2006) and earlier papers, the 
Ku-band NRCS over ocean at 10° incidence angle is 
—7.3 dB as based on TRMM satellite and airborne mea- 
surements. This value is also close to what is calculated 


from ocean-surface scatter models. On the basis of this 
well-accepted NRCS, we adjusted our internal cali- 
bration to about 1 dBZ lower to match this value. The 
HIWRAP-derived NRCS was obtained using surface 
returns over the Gulf of Mexico, where the ER-2 per- 
formed roll maneuvers. For Ka band, we followed an 
approach similar to that of Tanelli et al. (2006) in which 
the Ka- and Ku-band reflectivities are matched near 
cloud top where the scatterers are assumed to be Ray- 
leigh scattering at both frequencies. We adjusted the Ka- 
band reflectivity downward by —1.7 dBZ to match the 
Ku band. This provides a consistent calibration ap- 
proach that resolves any fluctuations in transmit power 
or receiver gain that occur during flight as a result of 
changes in temperature and altitude. The final calibra- 
tion should have an accuracy of better than 1 dB at Ku 
and Ka bands relative to TRMM and APR-2 and, it is 
presumed, close to an absolute calibration. 

The folding and accuracy of the Doppler velocities are 
of interest because the vertical motions in the sub- 
sequent data are very large. HIWRAP uses dual pulse 
repetition frequencies (PRF) to expand the Nyquist in- 
terval. For the 3589- and 4516-Hz PRFs used at both 
frequencies, the corresponding Doppler Nyquist veloc- 
ities are 20 and 25 ms -1 at Ku band and are 9 and 
12 ms -1 at Ka band (Table 1). With the staggered PRF 
approach used by HIWRAP, the effective Nyquist ve- 
locities are a larger by a factor of 4 than the single-PRF 
estimates above. These staggered PRF estimates were 
used to unfold the single-PRF Doppler velocity esti- 
mates since the latter have lower standard deviations 
(e.g., Doviak and Zrnic 1993, 171-175). This approach 
worked well for Ku- and Ka-band velocities. The accu- 
racy of the Doppler velocities is discussed in the ap- 
pendix as it relates to this paper. 
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Fig. 1 . Horizontal cross section of reflectivity at 3 km AGL from KOUN at 0121 UTC 24 May 
and KVNX at 2143 UTC 24 May. The red arrows are the ER-2 flight tracks during 0122-0126 
and 2145-2148 UTC. The locations of the OUN and LMN soundings and the KOUN and 
KVNX radars are marked. 


3. Overview of the hailstorms 

Conditions were highly favorable for severe weather 
during 23-24 May 2011 in the Oklahoma-Kansas region. 
There were two ER-2 flights and several severe thun- 
derstorm overpasses during this period; the first flight 
was from 2133 UTC 23 May to 0131 UTC 24 May, and 
the second flight was between 1920 and 2240 UTC on 
24 May. In the following, we will give a general description 
of the environmental conditions and time evolution of 
the storm. Two cases of intense convective storms are 
presented in this paper. Case 1 (Cl) during the first flight 
was located northwest of the S-band polarimetric ra- 
dar located at Norman, Oklahoma (KOUN). Case 2 
(C2) during the second flight was located near the 
Kansas-Oklahoma border within range of the Weather 
Surveillance Radar-1988 Doppler (WSR-88D) S-band 
polarimetric radar located at Vance Air Force Base in 
Oklahoma (KVNX). Figure 1 shows the general loca- 
tions of the two flight lines, S-band radars, soundings, 
and storms of interest superimposed on a 3-km AGL 
reflectivity map. All heights in this paper are above 
ground level, with the exception of upper-air sounding 
heights, which are above mean sea level. 


a. Upper-air soundings 

The upper-air soundings at Norman (OUN) and Lamont 
(LMN), Oklahoma, provide general environmental con- 
ditions for Cl and C2, respectively. The 0000 UTC OUN 
sounding (Fig. 2) was closest to Cl and showed a pseu- 
doadiabatic convective available potential energy 
(CAPE) of —3600 J kg -1 , but special soundings launched 
as part of MC3E at Purcell, Oklahoma, showed even 
larger instability, with CAPE of 3800-4200 J kg -1 be- 
tween 1732 and 2030 UTC. The 1800 UTC LMN sounding 
(Fig. 3) was closest to C2 and showed a CAPE of 
-4400 J kg -1 ; the 2100 UTC LMN sounding had a 
slightly lower CAPE of 3600 J kg - 1 . The freezing level is 
at —4. 3-km altitude in both soundings. The extreme 
instability in these soundings provides a ballpark esti- 
mate of over 80 m s~ 1 using iv max = (2 X CAPE) 172 , which 
indicates the maximum vertical motion at the equilibrium 
level calculated from the positive area in CAPE. The 
equilibrium level is — 175 hPa (— 13 km MSL) for the 
above soundings. The CAPE-derived updraft will never 
be realized because of precipitation loading and other 
factors, but it suggests the likelihood of very intense up- 
drafts. The bulk Richardson number was 41m 2 s -2 at 
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Fig. 2. Upper-air sounding for storm Cl from OUN at Fig. 3. As in Fig. 2, but from LMT at 2100 UTC 24 May 2011 for 
0000 UTC 24 May 2011. The dotted line shows the lifted-parcel storm C2. 

pseudoadiabat. 


OUN (Fig. 2) and 107-213 m 2 s -2 at the LMN soundings. 
These values are higher than typical supercell environ- 
ments (e.g., Weisman and Klemp 1982). 

b. Storm Cl 

During Cl, an intense convective line developed to 
the northwest of KOUN with anvils trailing eastward. 
Figure 4 shows a time evolution of the peak reflectivity 
in the Cl flight line as derived from KOUN plan position 
indicator (PPI) reflectivity scans. The maximum KOUN 
reflectivities exceed 70 dBZ prior to the ER-2 pass, with 
values as large as 75 dBZ around 15-20 min prior to the 
pass. The 60-dBZ echo extends up to 12-13 km throughout 
the storm’s period, and the 55-dBZ contour extends up to 
15-km altitude 5-10 min prior to the ER-2 pass, suggesting 
a strengthening of the updraft. 

The ER-2 made several passes over large cells during 
the evening of 23 May and into the early morning of 24 
May. The maximum reflectivity observed by KOUN 
is 70 dBZ for a cell located a few kilometers north of 
the flight line. By the time the ER-2 flew over Cl at 
0121 UTC, the reflectivity decreased slightly but still 
exceeded 65 dBZ. Figure 5 shows that the differential 
reflectivity Z DR at 3 km AGL within the high-reflectivity 
region was near or slightly below OdB. Such a feature, 
called the Z DR hail signature by Bringi et al. (1986) and 
Aydin et al. (1986), has been used to detect hail shafts 
penetrating below the melting level (Fig. 2). There also 
appears to be a “Z DR ring” signature that is frequently 


observed in supercells in midlevels (e.g., Kumjian and 
Ryzhkov 2008; Payne et al. 2010). The horizontal cross 
section at 7 km (Fig. 6a) shows a Z DR column located 
about 1 km south of the flight line reaching up to 7 km. 
This is accompanied by a reduction of cross-correlation 
coefficient p hv (Fig. 6b), which may be an indication of 
wet growth of large hail (diameter > 2.5 cm as defined by 
the U.S. National Weather Service) (Picca and Ryzhkov 
2012). There were numerous hail reports in the National 
Oceanic and Atmospheric Administration Storm Pre- 
diction Center storm reports, including five within 15 min 
of the radar data shown in Fig. 5 with hail sizes between 
2.54 and 5 cm. The ER-2 flight lines over this storm are 
about 80 km away from the KOUN radar so that its beam 
resolution is about 1.3 km wide with a 250-m gate spacing, 
which means that hail detection should not be a problem 
given the large size of the high-reflectivity core. 

c. Storm C2 

Later in the day on 24 May, storms first fired up at 
about 1840 UTC along the dryline, and then they de- 
veloped into a north-south-oriented line of supercells 
during the latter part of the second ER-2 flight. A pos- 
sible weak tornado, hail reports, and strong winds were 
associated with this group of storms near the time of the 
C2 overpass. After 2200 UTC, there was at least one 
report of 4-5-cm hail and several possible tornadoes 
with these same storms, as well as one report of 10-cm 
hail in the vicinity of C2. At least 12 tornadoes de- 
veloped with this line of supercells, mainly to the south. 
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Fig. 4. Time history of maximum reflectivity for storm Cl overflown by ER-2 (Fig. 1). The 
vertical dashed lines bound the period of the ER-2 flight over the storm. Contour levels are for 
50, 65, and 70dBZ (lighter to darker); calculated maximum values are also indicated (dBZ). 
The plot was derived from KOUN WSR-88D PPI scans. 


Figure 7 reveals that reflectivity values were not as 
high as for Cl. The reflectivity did not exceed 70 dBZ 
except at approximately 2200-2208 UTC at which time 
a >70-dBZ echo appears at 7-9-km altitude and then 
descends toward the surface afterward. The 50-dBZ 
contour extends up to about 11km as compared with 
13 km for Cl. The ER-2 flew across an intense convec- 
tive cell north of KVNX (Fig. 8). Unlike the previous 
case, we did not observe the Z DR ring and hail signature 
near the surface. Although a positive Z DR column 
reached up to 5 km AGL, we did not see an obvious 
reduction of p/„, for this particular cell (not shown). In 
the following section, we present the observations from 
Cl and C2 to gain an understanding of both the vertical 
storm structure and the response of the Ku and Ka bands 
to intense hailstorms. 

4. HIWRAP and ground-based radar observations 

of hailstorms 

a. Storm Cl flight line 
1) Vertical cross section along ER-2 

FLIGHT LINE 

To facilitate comparison between the HIWRAP Ku-/ 
Ka-band data and the S-band KOUN data, we have 
mapped KOUN data from scans in spherical radar co- 
ordinates to latitude, longitude coordinates for the ER-2 
flight line as previously used in Heymsfield et al. (2000) 
and Tian et al. (2002). Figure 9a shows a vertical cross 
section of S-band reflectivity from KOUN volume scans 
between 0121 and 0125 UTC along the ER-2 flight track 


during 0122:44-0126:19 UTC (shown in Fig. 5). The 
storm top exceeded 16 km in both the KOUN scan 
(Fig. 9a) and HIWRAP reflectivity (Figs. 9b, c). An in- 
tense core observed by KOUN with reflectivity greater 
than 60 dBZ and a horizontal dimension of about 10 km 
extends from the surface up to 11-km altitude. This core 
is sloped toward the southeast with height and is con- 
sistent with the northwesterly surface-500-hPa vertical 
shear (—25 ms -1 toward 130°) during the period. The 
motion of Cl was —4.4 ms -1 toward 150° so that most of 
the observed tilt was due to shear rather than from storm 
advection during the time taken for the PPI scan se- 
quence used in constructing Fig. 9a. This deep core of 
high reflectivity along with near-zero Z DR reaching the 
surface is a clear indication of a hail shaft (Fig. 9d) (e.g., 
Bringi et al. 1986; Wakimoto and Bringi 1988). A Z DR 
column up to 7.5 km to the south of the flight line sug- 
gests that liquid and mixed-phase hydrometeors were 
lofted to high altitudes, as will be discussed later. In this 
same region, HIWRAP reflectivity has significant at- 
tenuation. The reflectivities from the Ku and Ka bands 
are completely attenuated below 8 and 12 km, respec- 
tively (Figs. 9b, c). The peak reflectivity through this 
region is —45 dBZ for Ku band and —35 dBZ for Ka 
band as compared with the 63 dBZ observed at S band. 
This is not surprising since large hail would be in the 
Mie-scattering and highly attenuating region at Ku and 
Ka bands, as will be discussed further in section 5. 

The large upward Ku-band Doppler velocities from 
HIWRAP (Fig. 9e) have been corrected for aircraft 
motion. Doppler velocities by convention are positive 
downward (away from HIWRAP), but we have swapped 
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Fig. 5. Florizontal cross section at 3 km AGL at 0121 UTC 24 May. The black contours are 
KOUN reflectivity starting at 30 dBZ at 10-dBZ intervals. The Z DR values are color shaded, 
starting at — ldB. The red line is the ER-2 flight track between 0121 and 0126 UTC. The 
maximum reflectivity is 65 dBZ and is marked by a plus sign. Flight-line direction is from 
southeast to northwest. 


the sign on the color bar so that positive Doppler ve- 
locities are upward. The Ka-band Doppler velocity (Fig. 
9f) shows a structure that is similar to that of the Ku 
band in general except in the storm’s core. A strong 
upward Doppler velocity core above 10 km AGL at 
~7-km distance in Fig. 9e suggests an extremely strong 
updraft when the fall speed contribution is added (see 
section 5); Fig. 9f shows that the top of this updraft is 
above 10-km altitude before attenuation, Mie scattering, 
and multiple scattering become significant. The Doppler 
velocities in this updraft profile are not continuous from 
the maximum near 11 km to the surface, which may be 
due either to the tilt of the updraft in or out of the ver- 
tical plane of the measurements or to large hail fall 
speeds below the updraft maximum. At Ku band, the 
signal is strong in this region and therefore the Doppler 
estimate should be reliable (see the appendix); Ka-band 
Doppler has more problems in this region, with many 
factors contributing to poor measurements. This situation 
includes even the possibility that the Ka-band Doppler 
is affected by multiple scattering similar to what was 
modeled by Battaglia et al. (2011) for the W band. 
Strong downward hydrometeor motions (<— 30ms -1 ) 


at higher altitudes are present to the northwest of the 
updraft core (distance 11-16 km in Fig. 9e). The spatial 
resolution (—175 m at 10-km altitude) of the Ka-band 
Doppler measurements (Fig. 9f) is —2.5 times that of Ku 
band, and the Ka-band measurements show similar 
downward hydrometeor motions near this feature. 
Therefore, this feature is most likely due to a downdraft 
and/or high fall speeds rather than to other effects such 
as nonuniform beam filling (NUBF) that could bias the 
Doppler measurements (e.g., Tanelli et al. 2002). 

2) Averaged vertical profile 

We next examine the vertical structure near Cl’s core 
by averaging vertical profiles in the region between the 
two dashed lines in Fig. 9. Figure 10a shows mean pro- 
files of reflectivity and the difference of the measured 
reflectivity at two frequencies (reflectivity at lower fre- 
quency minus the reflectivity at higher frequency). The 
best agreement between the reflectivity at the three 
frequencies occurs near cloud top (>14 km) since the ice 
particles are small and closer to Rayleigh scatterers for 
all three frequencies, and attenuation is small. Note that, 
although reflectivity at S band is larger than 60 dBZ, the 
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Fig. 6. Horizontal cross section of (a) Z DR and (b) p hv at 7 km AGL at 0121 UTC 24 May. The black contours are KOUN reflectivity 

starting at 30 dBZ at 10-dBZ intervals. 


reflectivities at Ku and Ka bands at any level do not 
exceed —47 and —34 dBZ for this line, respectively. This 
difference is mostly due to Mie scattering at the very top 
of the high-reflectivity region and then a combination of 
Mie scattering, attenuation, and multiple scattering be- 
low. All of these effects are significantly worse at Ka 
band than at Ku band. As shown in Fig. 10a, the re- 
flectivities at Ku and Ka band are 35 and 60 dBZ, re- 
spectively, lower than that of S band at 6-km altitude. 
The Ka-band reflectivity becomes mostly attenuated 
below — 5-km altitude. Separating the effects of Mie 
scattering, attenuation, and multiple scattering in the 


current observations shown in Fig. 10a is difficult. 
Battaglia et al. (2011, their Fig. 4) show that multiple 
scattering artificially increases the reflectivity and could 
potentially enhance the value of the peak reflectivity. 
Battaglia et al. (2006) further suggest that multiple scat- 
tering could be an issue for GPM DPR at Ka band since 
the spaceborne radars have larger footprints. As a result, 
the enhancement could be as large as tens of decibels. 
HI WRAP has a much smaller footprint, and therefore 
the effect of multiple scattering is less significant. 

Figure 10b again shows the strong updraft core with 
a peak Doppler value of —39 ms -1 at 11-km altitude, 



Fig. 7. Similar to Fig. 4 except that the Dodge City, Kansas (KDDC), radar was used to ex- 
amine storm C2. 
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Fig. 8. Similar to Fig. 5, but at 2143 UTC 24 May and the black contours are KVNX re- 
flectivity and the red line is the ER-2 flight track between 2145 and 2148 UTC. The maximum 
reflectivity is 65 dBZ and is marked by an X. 


where the corresponding Ku-band reflectivity is —45 dBZ 
and the S-band reflectivity is —60 dBZ (cf. Figs. 10a and 
10b) in the updraft core region. 

3) Vertical cross section through Z dr 
column 

Figure 11 shows a vertical cross section of KOUN 
measurements reconstructed from KOUN PP1 volume 
scans between 0121 and 0125 UTC along the dashed line 
passing through the Z DR column located south of the 
ER-2 flight track shown in Fig. 5. Similar to Fig. 9a, the 
>60-dBZ reflectivity core extending to the surface as- 
sociated with near-zero Z DR indicates a hailshaft. The 
positive Z DR column at a distance of 6 km reaches up to 
7.5 km (Fig. lib) and is capped by a reduction of p hv 
(Fig. 11c). The higher value of Z DR of 4dB below the 
freezing level is probably due to the melting of hail- 
stones. As recent observational and modeling studies 
show, the Z DR column pinpoints the localization of 
a convective updraft, and its vertical extension is pro- 
portional to the updraft strength (e.g., Kumjian et al. 
2012; Picca et al. 2010). Supercooled raindrops are 
usually responsible for relatively shallow Z DR columns 
(1-2 km above freezing level, which is at —4.3 km in this 


case). Taller columns observed here usually signify that 
large graupel/hail particles are growing in the wet re- 
gime and that they are water coated, producing high 
Z DR at midlevels in the updrafts. Of course, such wet 
growth is not possible at very low temperatures, and 
upper parts of strong updrafts do not produce Z DR sig- 
natures at all. Wet growth usually manifests larger-size 
hail and is additionally marked by substantial depression 
of the cross-correlation coefficient at the levels at which 
the temperature is between —10° and — 20°C (Picca and 
Ryzhkov 2012). There is a strong indication that giant 
hail with sizes exceeding 5 cm is usually associated with 
strong p lw depression and slightly negative Z DR above 
the melting layer (Kumjian et al. 2012). 

b. Storm C2 flight line 

Figure 12 shows the vertical cross section along the 
flight line in Fig. 8. The radar cloud top for this storm 
extended up to 15 km AGL, which is —2 km lower than 
for Cl. Figure 12 indicates the appearance of a weaker 
storm than Cl. The maximum reflectivities are high but 
barely reaching 60 dBZ at S band (Fig. 12a). The Ku- 
and Ka-band reflectivities (Figs. 12b, c) show an in- 
teresting “cap”-like reflectivity structure at —12 km 


1860 


JOURNAL OF APPLIED METEOROLOGY AND CLIMATOLOGY 


Volume 52 



0 10 20 30 40 

DIST(km) 


60 

50 


40 N 
CD 

302. 


U) 

20 N 

10 

0 

60 

50 

40 n 
CD 
30 2. 

20 n 


10 

0 

60 

50 

40 n 
CD 

30 B 
ro 

20 m 


10 

0 





15 

E 


h- 

I 

O 

10 

ID 

I 

5 


40 




Ku, HIWRAP 


D 40- 

Ka, HIWRAP 


0 10 20 30 40 

DISTIkml 



1 ° 
' N 

0 


-1 

2 


-40 


-30 


-20 


-10 

E, 

0 

Q> 

10 

Q. 

CL 

20 

a 

30 


40 



-40 


-30 


-20 

(/) 

-10 


0 

q3 

10 

Q. 

CL 

20 

a 

30 


40 



Fig. 9. Vertical cross section along the flight line at 0121-0126 UTC 24 May 2011, as shown in Fig. 1, for (a) S-band reflectivity from 
KOUN, HIWRAP reflectivities at (b) Ku and (c) Ka bands, (d) Z DR with KOUN reflectivity contoured, and HIWRAP Doppler velocity 
at (e) Ku and (f) Ka bands. The contour lines are S-band reflectivity starting at 30dBZ with an interval of 10 dB. The HIWRAP ob- 
servations extend up to approximately 16-km altitude because of the pulse compression used. The section is oriented from southeast to 
northwest. The plus sign in (d) indicates the location of the maximum reflectivity. The dashed vertical lines in (a)-(c) indicate the 
averaging region used in Fig. 10, below. The dotted vertical line in (b) is the center of this region. 


AGL. The reflectivities at Ku and Ka bands are com- 
pletely attenuated below 3 and 5 km AGL, respectively. 
The Z DR does not have a distinctive hail signature 
(Fig. 12c) at lower levels during the overpass time. A 
strong positive Doppler velocity is noted in the storm’s 
core that is associated with the updraft (Figs. 12, f), al- 
though these Doppler velocities are weaker than those 
in Fig. 9. The profiles in Fig. 13 also indicate that this 
storm is weaker than the Cl storm. The S-band re- 
flectivity (Fig. 13a) barely reaches 55 dBZ, the Ku-band 
reflectivity peaks at —47 dBZ, and the Ka band peaks at 
—37 dBZ at 12-km altitude. The lower S-band reflec- 
tivity relative to Ku above 1 1 km is mainly due to the 
differences in volume-scan and overpass times and 
spatial resolutions between HIWRAP and KVNX. 
Doppler velocities are also somewhat different from Cl, 
where peak upward motions (without fall speed added) 
of 25 m s” 1 occur at a lower altitude between 4 and 8 km. 

Several regions of downward Doppler velocities 
(<0 m s” 1 ), indicative of downdrafts or large fall speeds, 
are present at higher altitudes surrounding the main 


updraft region of C2. These were higher than for Cl, 
with a minimum just lower than —30 ms -1 at —10.5 km 
AGL and 22-km distance. Similar to what was observed 
in Cl, the Ka-band Doppler (Fig. 12f) showed down- 
draft features that are similar to those of the Ku-band 
Doppler (Fig. 12e), and therefore NUBF Doppler biases 
do not appear to be responsible for this feature. 

5. Interpretations from the observations 

a. Summary of observations 

The analysis of MC3E hailstorms with HIWRAP and 
polarimetric measurements has shown several prom- 
inent features: 1) strong updrafts are suggested by the 
Doppler velocities, with 39 ms -1 at 10.7km (250 hPa) 
and 30ms” 1 at 13.7km (160 hPa) in storm Cl, 2) large 
attenuation is observed in Cl and C2, 3) a tall Z DR 
column up to 7.5 km is accompanied by a depression of 
p hv , 4) S-band reflectivities exceed 70 dBZ, with 60 dBZ 
extending up to 10.5-km altitude, and 5) the Ku- and 
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Fig. 10. Profiles corresponding to Fig. 9: (a) averaged profile of 
reflectivity at three frequencies (S, Ku, and Ka) and reflectivity dif- 
ferences between the two frequencies (S — Ku, S - Ka, and Ku - Ka) 
and (b) Ku-band Doppler velocity at ~5-km distance in Fig. 9, along 
with the minimum and maximum Doppler velocities within the same 
section. Average profiles use 30 profiles in the section between the 
two vertical dashes in Fig. 9. The Ku-band reflectivity at the location 
of maximum Doppler velocity is also plotted in (b) (blue line). 

Ka-band reflectivities in Cl and C2 decrease signifi- 
cantly below 12-14-km altitude, relative to the S-band 
reflectivity. There are a number of questions related to 
both the storm microphysics and the use of spaceborne 
radar. Here, we address two of them: 1) What are the 
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DIST (km) 


Fig. 11. Vertical cross sections of (a) S-band reflectivity, (b) Z DR , 
and (c) p lw from KOUN volume scans between 0121 and 
0125 UTC along a line 1.5 km south of the ER-2 flight line. For 
reference, horizontal dashed lines show 7-, 8-, and 9-km heights. 
The black contours show the 60-dBZ reflectivity. 


magnitudes of the updrafts? 2) Are the ice particles at 
the higher altitudes graupel, or are they small or large 
hail? In the following, we discuss the current observa- 
tions to provide insight on these questions. 

b. Updraft magnitudes from Doppler observations 

Reflectivity-weighted fall speeds are required to esti- 
mate the vertical motions from the Doppler measure- 
ments. Ulbrich (1977) calculated reflectivity-weighted 
hail fall speeds for a range of wavelengths from S band 
(10 cm) down to X band (3.2 cm), for dry and wet hail 
with density 0.917 gem -3 , and for representative expo- 
nential size distributions that were based on empirical 
observations. A simplified v T = A'Z b e is given for dry or 
water-coated hail by Eq. (11) in Ulbrich (1977), where 
A' and b are coefficients in his Table 2, and Z e is the 
equivalent reflectivity. Coefficient A' is normalized to 
N 0 = 1 00 m in the exponential distribution so that it is 
not sensitive to its fluctuations. Ulbrich states that the v T 
is not very sensitive to wavelength for dry and water- 
coated hail. For the 60-dBZ S-band reflectivities near 
the Doppler velocity maximum in Fig. 9a, S-band v T for 
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Fig. 12. As in Fig. 9, but for a flight line between 2145 and 2148 UTC 24 May, as shown in Fig. 8. 


dry hail is 20.1ms -1 at the 1000-hPa level; the corre- 
sponding X-band v T is 20.9 ms -1 . Above the surface, v T 
must be corrected for decreased air density with altitude 
by an approximate factor of [p o lp(z)] 0A , where p(z ) and 
p 0 are the densities at height z and at 1000 hPa, re- 
spectively (Beard 1985). This would suggest an increase 
in Ulbrich’s fall speed estimate by a factor of —1.5 or 
—30 ms -1 at 250 hPa (10.7-km altitude). The above 
correction was based on rain, and more recent work has 
shown an exponent closer to 0.4 for snow and graupel so 
that there might be a small reduction in fall speed 
magnitude (A. Heymsfield 2012, personal communica- 
tion). A fall speed of 30 ms -1 and maximum Doppler 
velocity of 39 ms -1 would give a maximum updraft of 
69 ms -1 (Doppler velocity + v T ). 

Particles with lower density will have a smaller fall 
speed when compared with solid ice. The World Me- 
teorological Organization defines hail as solid ice 
(0.91 gem -3 ) and generally larger than 5mm, whereas 
graupel has a lower density (<0.4 gem -3 ) and is less 
than 5 mm; graupel is sometimes referred to as small 
hail. There are a number of reports in the literature of 
lower-density hail with values as low as 0.4 gem -3 for 
1-cm hail (e.g., Knight and Heymsfield 1983). If lower- 
density hail were present, it would still have a significant 


fall speed (e.g., Fig. 6 in Knight and Heymsfield 1983), 
although it is complicated to assess the overall impact on 
v T given that the size distribution is unknown. Reducing 
the v T estimate by half would still provide an updraft of 
54ms -1 . 

There are limited data to corroborate these peak up- 
draft estimates. Musil et al. (1991) presented measure- 
ments at 6-km altitude from an armored T-28 aircraft in 
severe Montana hailstorms, and they estimated up to 
53 ms -1 updrafts and 18 ms -1 downdrafts (Table 2 in 
Musil et al. 1991). Their updraft and downdraft core 
sizes tended to increase with updraft magnitude; the 
53 ms -1 updraft case had a core 15 km wide that con- 
tained 5-cm hail. Their observed updrafts tended to be 
stronger than downdrafts. Their most intense case on 
2 August 1981 had 5-cm hail with a hailstone concen- 
tration of 30 hailstones per meter cubed. Bluestein et al. 
(1988) studied a dryline supercell in the Texas Panhandle 
that was tornadic and producing hail. The updraft speeds 
estimated from sounding ascent rates were 35-40 ms -1 
between 6- and 7-km altitude, with a peak value of 
38ms -1 at 8km. Nelson (1983) calculated a 50ms -1 
updraft in an Oklahoma hailstorm from a dual-Doppler 
analysis. The updraft magnitudes (54-69 ms -1 ) obtained 
in this study are also consistent with the large CAPE in 
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Fig. 13. As in Fig. 10, but for a line between 2145 and 2148 UTC 
24 May. 


the environment (section 3) and the nearly 4-km over- 
shoot of the equilibrium level by the storm top. 

c. Hail scattering and attenuation 

Airborne radar measurements at Ku and Ka bands 
from APR-2 have been used for testing spaceborne 
dual-wavelength rain-retrieval algorithms (e.g., Liao and 
Meneghini 2011; Grecu et al. 2011; Le and Chandrasekar 
2012). Most of these studies have focused on stratiform 


rain because of its relatively simple vertical structures. 
The scattering properties of snow aggregates above the 
melting layer, mixed phase inside the melting layer, and 
raindrops below are better understood from field ex- 
periments and theoretical computations. Comparable 
studies have not been made for hail at Ku and Ka bands, 
mainly because of the lack of observations in severe 
storms, especially over land with good coverage of 
ground-based radar. In the following, we provide some 
insight on the hail reflectivities for HI WRAP and GPM 
and use the WSR-88D S-band wavelength as a refer- 
ence, similar to the presentation in section 4. We are 
particularly interested in whether the Ku- and Ka-band 
observations can be used to detect the presence of hail in 
future satellite measurements. 

The focus here is on obtaining a general sense of the 
behavior of large hail in the 10-12-km-altitude region of 
the storms where the HIWRAP reflectivities at Ku and 
Ka bands undergo rapid departures from the S-band 
reflectivities and where temperatures are sufficiently 
low. To understand fully the vertical profiles at Ku and 
Ka bands in such a hailstorm, we need to consider not 
just hail but also graupel with different shapes, densities, 
and particle size distributions, which is beyond the scope 
of this paper. The assumption is made that dry hail with 
density 0.917 g cm” 3 is present in this region. 2 Scattering 
calculations are performed at the three frequencies (S, 
Ku, and Ka band) using Mie theory with the assumption 
of a monodisperse distribution of dry spherical hail- 
stones up to 6 cm and a temperature of — 30°C. Figure 14 
shows calculations of the radar reflectivity, the one-way 
attenuation coefficient normalized by number concentra- 
tion N, and the dual-frequency ratio (DFR = Zfl ~ Z f2, 
where the subscripts indicate the radar frequency, with 
/ 1 < / 2). DFR is independent of number concentration. 
Given that we do not know the hail sizes or concentra- 
tions within the updraft regions, Fig. 14 conveys the 
essential characteristics of the hailstone scattering. 

Figure 14a shows that reflectivities at the three fre- 
quencies depart significantly from each other for hail of 
>0.5 cm in diameter. The reflectivity of Ka band is lower 
than those at Ku and S bands. This is apparent in the 
observations, for which reflectivities are never higher 
than 35 dBZ at Ka band and 55 dBZ at Ku band. Note 
that the S-band reflectivities also depart from Rayleigh 


2 The environmental air temperature is approximately -40°C 
(— 45°C) from a representative sounding at 10-km altitude from 
OUN (LMN) at 0000 UTC (1500 UTC) on 24 May 2011 for the Cl 
(C2) overpass, and the undiluted lifted parcel temperature is 

- — 35°C. The expected air temperatures are sufficiently low to 
ensure that all particles are frozen above 10-km altitude. Most 
raindrops freeze at much lower altitudes. 
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Fig. 14. Mie calculation for monodispersed hail size distribution, 
assuming ice density of 0.917 gem -3 , for (a) ZfIN , (b) -Z/i — Z/2> 
and (c) kf/N as a function of ice particle size, where Zy, kf , and 
Zfi — Zf 2 are the reflectivity, one-way attenuation coefficient, and 
difference of the reflectivity at two different radar frequencies (f 1 > 
/2), respectively, and N is the concentration of monodispersed 
particles. In (a) and (c), solid black, red, and green curves are for 
Ku, Ka, and S bands, respectively. In (b), red, black, and green 
curves are for Z s — Z Ku , Z s — Z Ka , and Z Ku — Z Ka , respectively. 

scattering for hailstones of 2 cm and larger. The DFR 
(Fig. 14b) is more complex. For hail diameters larger 
than 2 cm, the DFR is larger than 10 dB for Z s — Z Ku 
and Z Ku — Z Ka , and 15 dB for Z s — Z Ka . Between 0.5 
and 2 cm, DFR has a negative minimum near 1 cm for 


Z Ku — Z Ka and a large 17-dB peak near 1-cm size. This 
result illustrates why size distribution knowledge is so 
critical for understanding the hail DFR. 

Attenuation at Ku and Ka bands (Fig. 14c) shows 
a much more complicated behavior, with the two curves 
crossing at about 1, 3, and 4 cm. For hail sizes of less than 
~1 cm, attenuation at Ka band is larger than at Ku band. 
This trend is reversed for hail sizes between —1.2 and 
3 cm, however, and between 4.5 and 5.5 cm. Note that for 
hail sizes of —4 cm attenuation between the Ku and Ka 
bands are comparable. 

Large dry hailstones with diameters of more than 5 cm 
attenuate more at S band than at Ku and Ka bands. This 
is because at these sizes attenuation is caused mostly by 
scattering rather than by absorption. Attenuation at S 
band in hailstorms is rarely observed (A. Ryzhkov 2012, 
personal communication), however, because of the sig- 
nificantly higher concentration of smaller hailstones or 
graupel for which attenuation at S band is much lower. 
Therefore, the integrated attenuation at S band is dra- 
matically smaller. In a similar way, hailstones with di- 
ameters larger than 3 cm attenuate more at Ka band 
than at Ku band, but observed Ka-band attenuation is 
higher than that at Ku band. Assuming different hail 
densities or water coating on the hail will greatly modify 
the behavior of these curves. 

Figure 15 shows the Z Ku — Z Ka for the Cl and C2 
overpasses, where regions with >14dB are shaded 
white. The most striking feature is that the large DFR 
values are near 14 km (12 km) for Cl (C2). The highest 
portions of this high-DFR region (>14dB) are in the 
regions of strongest inferred updrafts from the Doppler 
velocities. The updraft regions likely have the highest total 
ice mass. The extent to which the high-DFR region in- 
dicates hail size and whether graupel is present is ambig- 
uous, however, because of the complex nature of the hail 
scattering and attenuation, as mentioned previously. It 
would take extremely high concentrations of graupel to 
account for the 50+ dBZ observed above 10-km altitude at 
Ku band in the storm cores. It is certainly possible, 
though, that there are mixtures of hail and graupel in this 
region. What is also of interest is that these large DFR 
values also extend to lower altitudes. For example, the 
cell in Fig. 9b (distance 25 km) has a large DFR value 
present up to 8-km altitude. In those regions, super- 
cooled cloud drops, unfrozen rain, graupel, or recircu- 
lated hail could be present at this level (— 20°C). Certainly 
any mixed phase, wet aggregates, or graupel at lower 
altitudes would cause a large DFR. Differential atten- 
uation between the two wavelengths contributes to this 
persistent large DFR, especially farther into the core. 
Further study is needed to determine the usefulness of 
DFR at Ku and Ka bands for hail detection. 
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Fig. 15. The Ku-/Ka-band DFR (Z Ku - Z Ka ) for (top) Cl and 
(bottom) C2. The white-shaded region shows DFR >14 dB. 


6. Conclusions 

New dual-frequency (Ku and Ka bands) nadir-pointing 
radar measurements from HIWRAP mounted on the 
ER-2 aircraft have been obtained from severe hail- 
storms in Oklahoma and Kansas. These measurements 
are at frequencies that are similar to those of the TRMM 
PR (Ku band) and the upcoming GPM DPR (Ku and Ka 
bands). For GPM measurements, it is important to un- 
derstand the dual-frequency radar characteristics in 
a wide variety of thunderstorm types and intensities 
globally. There is a wide range of vertical structures in 
convection because of varying updraft characteristics, 
microphysical composition such as hail or graupel, and 
their horizontal dimensions. This paper looks at more- 
extreme convection cases in the southern Great Plains. 
What are the limitations of spaceborne-radar measure- 
ments from GPM in strong midlatitude convection and 
what information can we expect to retrieve? Aircraft 
overflights of two strong hailstorms have provided not 
only storm structure but also implications of these kinds 
of storms for satellite measurements. Updrafts in severe 
storms have been previously estimated from S-band 
ground-based multi-Doppler analyses, inferred from ob- 
served soundings, or — less frequent — obtained from di- 
rect in situ aircraft measurements. Data from storms Cl 
and C2 have inferred updrafts possibly exceeding 
50-60 ms -1 for Cl and >35 ms -1 for C2. The Cl updraft 
magnitude is slightly higher than the strongest updrafts 
reported in the literature. The observed updraft is capable 
of supporting multicentimeter-size hail at higher altitudes. 
The Z DR column height (7.5 km AGL, or —3.2 km above 


the freezing level) corresponding to extreme updraft 
speeds is consistent with the theoretical modeling of 
Kumjian et al. (2012). 

The reflectivities at Ku and Ka bands from HIWRAP 
exhibited large Mie scattering, high attenuation rates, and 
likely multiple scattering in cores of both storms. As a 
result, the maximum reflectivities observed were 35 dBZ 
at Ka band and 45 dBZ at Ku band, whereas the values 
at S band were sometimes as high as 70 dBZ in the same 
regions. This is important for TRMM or GPM mea- 
surements since the peak heights of reflectivity levels 
such as 50 dBZ are often used to infer updraft strength in 
the absence of Doppler measurements. This procedure 
may result in ambiguity in interpretation since a 50-dBZ 
echo height at Ku band may actually be a much higher 
reflectivity. Calculations with simplistic hail assumptions 
(dry, spherical, and monodisperse) show that Mie effects 
become more significant for hail sizes of approximately 
1 cm and above for Ku- and Ka-band measurements. The 
DFR (Ku — Ka) exceeds 14 dB near the top of the high- 
reflectivity region at S band (10-12-km altitude) where 
Mie effects, attenuation, and multiple scattering all are 
likely significant. While this DFR information in itself 
may not indicate a unique determination of the ice 
particle characteristics (hail or graupel), it may be useful 
for the spaceborne measurements for differentiating 
a 50-dBZ storm from a 70-dBZ storm. We have not fo- 
cused on DFR interpretation at lower altitudes in the ice 
region where there are large variations in hydrometeor 
type (snow, graupel, hail, cloud, and supercooled rain- 
drops) because of their complexity. 

The horizontal extent of the SGP hailstorms is often 
10-20 km across and would be detected by the GPM 
DPR with the 5-km GPM DPR spatial resolution. If 
processes such as multiple scattering are indeed impor- 
tant in the aircraft measurements, they will be much 
worse in the measurements from GPM because of its 
larger footprints. Further modeling calculations will be 
required to assess these kinds of complications, as well 
as the general utility of dual-frequency measurements in 
convection. Future work will examine a variety of con- 
vective intensities along with more-realistic modeling of 
the scattering properties such as in the modeling 
framework provided by Tanelli et al. (2011) that deals 
with single and multiple scattering. 
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APPENDIX 

Accuracy of Doppler Measurements in Updraft 
Regions 

For Doppler measurements, adjacent pulses must be 
correlated. The decorrelation time is given by t, = 2A la v 
[Sauvageot 1992, his Eq. (1.104)] and must be larger 
than the interpulse time T s = 1/PRF, where A is the radar 
frequency, cr v is the target spectral width, and PRF is the 
pulse repetition frequency. If we use cr v = 10 ms -1 as 
our worst case, then r, is —4.4 ms at Ku band and 1.8 ms 
at Ka band. Since T s is 0.22 ms for both HIWRAP fre- 
quencies at the higher PRF (4516 Hz), the adjacent 
pulses are correlated as required for a Doppler mea- 
surement at Ku and Ka bands for the more severe 
conditions in the updraft cores. Note there is also de- 
correlation that is due to the aircraft motion, expressed 
by Tj = 0.57A/( Wl/j) [Meneghini and Kozu 1990; their 
Eq. (2.82)], where V is the platform motion and is the 
half-power antenna beamwidth from Table 1. This 
decorrelation time is —1.4 ms (1.8 ms) for Ku band (Ka 
band), assuming an aircraft motion of 160ms -1 so that 
adjacent pulses are still correlated at both frequencies. 

Evaluation of cr v from the observations is required to 
understand the magnitude of the errors in Doppler 
measurements. Various factors such as turbulence, shear, 
platform motion, and particle fall speeds all contribute to 
cr v , especially in the most turbulent regions such as in 
intense updrafts. The a v is calculated according to Eq. 
(6.27) in Doviak and Zrnic (1993), using estimates of the 
autocorrelation coefficient and the signal after noise 
removal: 
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where S is the signal power and R\ is the covariance at lag 1. 


Melnikov and Zrnic (2004) found that this 0,1-lag 
autocovariance estimator does better than other esti- 
mators for the large spectral widths in which we are in- 
terested here. They suggest the maximum spectrum 
width must be less than about D.6v N , where v N is the 
Nyquist velocity. This maximum width can also be ex- 
pressed by (v n /tt) ln(M 0 ' 5 ), where M is the number of 
samples. For the HIWRAP observations presented in 
this paper, u^is —24 ms -1 for Ku band and 10 ms -1 for 
Ka band, assuming a maximum PRF of 4516 Hz and 64 
samples. This implies that the spectral width estimator is 
valid for cr v < 15.9 ms -1 for Ku band and a v < 6.6 ms -1 
for Ka band. In the more turbulent updraft regions of 
the flight lines presented in this paper, spectral width 
estimates from Eq. (Al) peaked at —10 ms -1 in the Ku- 
band data. The Ka-band spectrum widths were un- 
realistically low in the updraft region, suggesting that the 
Eq. (Al) assumptions were failing (figure not shown). In 
general, we would expect lower values at Ka band 
because of the narrower beamwidth of the HIWRAP 
Ka-band antenna as compared with Ku band. Our 
conclusion from the above is that the Ku-band ve- 
locities are valid in the convective updrafts since ad- 
jacent pulses are correlated based on spectral width 
estimates from the data, and the spectral widths are 
less than QAv N even in the updraft regions. The large 
spectral widths combined with low signal-to-noise 
ratio in highly attenuated updraft regions are likely the 
cause of poor or missing velocity estimates at Ka band, 
however; this is the case in Figs. 9f and 12f. The variance 
of the Doppler velocity can be estimated through per- 
turbation analysis. For large signal-to-noise ratios 
and a narrow cr v = 10 ms -1 , var(u) — crA/(8Mr s 7r 1/2 ) 
[Doviak and Zrnic 1993, their Eq. (6.23)] provides an 
uncertainty of — 1 m 2 s -2 for Ku-band Doppler velocity. 
We have not provided corresponding Ka-band error 
estimates in these worst-case situations since the autoco- 
variance estimators do not perform well at this frequency. 
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